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Abstract 

The form factors of the semileptonic Bg — > Dq{J^ = 0^)£i' with 
q = s,u,d transitions are calculated in the framework of three point 
QCD sum rules. Using the dependencies of the relevant form fac- 
tors, the total decay width and the branching ratio for these decays 
are also evaluated. A comparison of our results for the form factors 
oi B ^ Diiy with the lattice QCD predictions within heavy quark 
effective theory and zero recoil limit is presented. Our results of the 
branching ratio are in good agreement with the constituent quark me- 
son model for (g = s,u,d) and experiment for (g = u,d). The result 
of branching ratio for Bg Ds(1968)£i^ indicates that this transition 
can also be detected at LHC in the near future. 
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1 Introduction 



The pseudoscalar Bq meson decays are very promising tools to constrain 
the Standard Model (SM) parameters, explore heavy quark dynamics and 
search for new physics. The semileptonic decays of heavy flavored mesons 
are also useful for determination of the elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, leptonic decay constants as well as the origin of the 
CP violation. Neutral and B^ meson decays are interesting to study CP 
violation. 

When LHC begins operation, an abundant number of Bq mesons will be 
produced. This will provide a real possibility for studying the properties of 
the Bq mesons and their various decay channels. Some possible decay chan- 
nels of Bq mesons are their semileptonic decays to Dqlu . The Bq Dqiu 
transitions occur via the b ^ c transition with s, d or u as spectator quarks. 
The most common decay mode of B mesons is clearly b ^ c transition, 
since it is the most dominant transition among the b quark decays. The 
semileptonic Bq — > Dqiu decays are interesting because they could play a 
fundamental role in probing new physics charged Higgs contributions in low 
energy observables. Moreover, they open a window onto the strong inter- 
actions of the constituent quarks of the pseudoscalar Ds meson and could 
give useful information about the structure of this meson (for a discussion 
about the nature of Dgj mesons and their quark content see [l, 2j). Analysis 
of the D.o(2317) ^ D:7, D.j(2460) -> D^-f and I).j(2460) ^ (2317)7 
indicates that the quark content of these mesons is probably cs p] . 

The long distance dynamics of such type transitions are parameterized in 
terms of some form factors, which are related to the structure of the initial 
and final meson states. For calculation of these form factors which play 
fundamental role in the analysis of these transitions, some nonperturbative 
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approaches are needed. Among the existing nonperturbative methods, QCD 
sum rules has received especial attention, because this approach is based 
on the fundamental QCD Lagrangian. There are two kinds of QCD sum 
rule approaches, three point and light cone QCD. In three point QCD sum 
rules, the perturbative part of the correlation function is expanded in terms 
of operators having different mass dimensions with the help of the operator 
product expansion (OPE). In light cone QCD, the distribution amplitudes 
(DA's) of the particles expanding in terms of different twists are used [H O 
[6]. This method has been applied successfully for wide variety of problems 
[3 El El [ini E] (for a review see also [12] )• In present work, we describe the 
semileptonic Bq Dgiu decays by calculating the relevant form factors in the 
framework of the three point QCD sum rules approach. Note that, the form 
factors of -B — s> Diu have been calculated in lattice QCD [T^l El [13, [IS] and 
the subleading Isgur-Wise form factor is computed in QCD sum rules and its 
application for the B Diu decay is shown in ^TTl [18] (for similar previous 
works see also [121 [201 [2T] )• Moreover, the Bq — >■ Dqii/ transitions have been 
studied in the constituent quark meson (CQM) model for q = s,u,d in [22] 
and for q = u,d, the experimental results can be found in [23]. 

This paper is organized as follows: In section II, we calculate the sum 
rules for the two form factors relevant to these transitions. Section III is 
devoted with the numerical analysis, conclusion, discussion and comparison 
of our results for the form factors and branching ratios with those of the 
other phenomenological model, lattice QCD and experiment. 
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Figure 1: Bq ^ Dqiu transitions at tree level 

2 Sum rules for the Bq Dqiu transition form 
factors 

In the quark level, the Bg Dgiu transitions proceed by the b ^ c transition 
(see Fig. 1). The matrix element for these transitions at the quark level can 
be written as: 

Mq = ^Kfe 7/.(l - 75)/ c 7^(1 - 75)&. (1) 

To obtain the matrix elements for Bq Dqiu decays, we need to sandwich 
Eq. ([1]) between initial and final meson states, so the amplitude of these 
decays gets the following form: 

M = ^V,t V 7^(1 - 75)/ < Dqijl) I c 7^(1 - 75)& I Bqip) > • (2) 

Our aim is to calculate the matrix elements < Dq{p') \ C7^(l —75)6 | Bq{p) > 
appearing in Eq. ([2]). Because of parity and Lorentz invariance the axial vec- 
tor part of transition current, c 7^(1 — 75)6 , does not have any contribution 
to the matrix element considered above, so the contribution comes only from 
the vector part of the transition current. Considering the parity and Lorentz 
invariances, one can parameterize this matrix element in terms of the form 



3 



factors in the following way: 

< D,ip') I C7,6 I B^{p) >= h{q^)P^ + Uq^)q,. (3) 
where f2{.<f') are the transition form factors and = {p + p')^, 

From the general philosophy of QCD sum rules method, in order to cal- 
culate the form factors we consider the following correlator: 

Mp^p", =^'J d'xdSe-'^^e'^'y < | T[Jz,,(i/) 7^(0) Jb,(x)] | >, (4) 

where iy) = c'j^q and Jb^ (x) = b'y^q are the interpolating currents of the 
Dg and Bg, respectively and ^7^(0) = 07^6 is the transition current. 

To calculate the phenomenological or physical part of the correlator given 
in Eq. (jlj), two complete sets of intermediate states with the same quantum 
numbers as the currents Jd, and Jb^ respectively are inserted. As a result of 
this procedure, we get the following representation of the above-mentioned 
correlator: 

< I Jp^jO) I Dgjp') >< Dgjp') I J^(0) I Bgjp) X Bgjp) \ JgjO) \ > 

{p'^-ml^){p'^-mlj 
+ (5) 



n 



where ■ ■ • represents the contributions coming from higher states and contin- 
uum. The following matrix elements in Eq. ([5]) are defined in the standard 
way as: 



<0\Jd,\ Dg{p') >-- 



. fD.ml 

— 

nic + ruq 



< Bg{p) I Jb, I >= -z (6) 

nib + nig 

where and fs^ are the leptonic decay constants of Dg and Bg mesons, 
respectively. Using Eq. and Eq. ([H]), Eq. (jS]) can be written in hadronic 
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language as: 
Where, 

nip,p,g) = -7^^ 7^7^ — \ — \ — -fm) 

f^piz _ rrij^^ ) [p"^ — rrijg^ ) nic + niq rrib + rriq 
+ excited states, 

^^2yp ,p ,q ) — —TI2 2 \r 2 2~T \ \ ^21,? ) 

[p'-^ — mo^yp — + + f^q 

+ excited states. (8) 

Now, let calculate the theoretical part (QCD side) of the correlation func- 
tion Iln{p^,p'^, q^) in quark and gluon languages with the help of the operator 
product expansion(OPE) in the deep Euclidean region p^ <^ {uib + rnqY and 
p'^ <^ {rric + mq)^. The correlator is written in terms of the perturbative and 
nonperturbative parts as: 



+ 



^^(p^p'^g^)+^r-^^'•(p^p'^g2) 

Urip\p'\q')+Ur-''^{p',p",q') 



%■ (9) 



To obtain the sum rules for the form factors, the two different representa- 
tions of n^(p^,p'^, q^) are equated. The theoretical part of the correlator is 
calculated by means of OPE, and up to operators having dimension d = 5, it 

is determined by the bare-loop (Fig. 2a) and the power correction diagrams 
from the operators with d = 3, < gg >, c? = 4, < qq >, d = 5, rn^ <qq> 
(Fig. 2b, 2c, 2d). In calculating the bare-loop contribution, we first write the 
double dispersion representation for the coefficients of corresponding Lorentz 
structures appearing in the correlation function as: 

■Qper _ _ f [ ds'- — , ' ^ -l- subtraction terms. (10) 

' (27r)2 i J {s - p^){s' - p'^) ^ ' 

5 




Figure 2: Diagrams for bare- loop and power corrections (light quark conden- 
sates) 

The spectral densities pi{s,s',q^) can be calculated from the usual Feyn- 
man integral (bare loop diagram in Fig. 2a) with the help of Cutkosky 
rules, i.e., by replacing the quark propagators with Dirac delta functions: 
p2^^2 —2-K5ijP' — w?), which implies that all quarks are real. After some 
straightforward calculations for the spectral densities corresponding to 
and we obtain: 

pi(s, s',g^) = NJo{s,s',q^)[2mbmq + 2mcmg-4:ml 

- 2A{A + ml- ml) - {A + B)u' - 2B{A' + ml- m^)], 

P2(s, s',g^) = NJo{s,s',q'^)[-2mbmg + 2mcmg + 2A{A + ml - rrig) 

+ {B-A)u'-2B{A' + ml-ml)], (11) 

where 

Io{s,s',q^) — -rrT-pTf ; 57; 

4A^/^(s, s , q'^) 
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A' 
A 

A 

B 

u 
u' 



+ + - 2sg2 - 2s' - 2ss', 
is' -ml + ml), 



[s-ml + ml), 



1 

s + s' -q"^ 



— -[25'A-A'«] 



[2sA'- Au], 



2[mbmc — {nih + mc)mq + m^]. 



(12) 



In Eq. ffTTl) A'^c = 3 is the number of colors. The integration region for the 
perturbative contribution in Eq. ( |TOl) is determined from the condition that 
arguments of the three 5 functions must vanish simultaneously. The physical 
region in s and s' plane is described by the following inequalities: 



^ ^ 2ss' + (s + s' — q^){ml — s — m^) + (m^ — m^)2s ^ 



(13) 



A^/2(m^, s, m^)A^/2(s, s', g^) 
From the above equation, it is easy to calculate the lower bound of integration 
over s' as a function of s. (i.e., s' = /(s)). 

For the contribution of power corrections, i.e., the contributions of opera- 
tors with dimensions d = 3, A and 5 (diagrams in Fig. 2b, 2c ,2d), we obtain 
the following results: 



non—per 



1 Tfi 

-—— <qq> {rrib + m^) + ^ <qq> 
Ir'r 4 

1 , o 1 



mbrric + ml ^ ml + m^mc 



+ ^<qq> {-ml + -ml 



mbml + ml 1 [m^ + mb){—q^ + ml + m\ 
^ 2 



1. 

-(mc + mfej 



1 1 

+ 



+ 



ml + mime 



ry ^2 fy 2 
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+ 



ml {Snib + rric) _ "^o / ~4(mb + 2mc) {2mb + rric) 
{—ml + mime + m^ml — m^ + m;,g^ + mcQ^) \ 



1 T/1 
jjnon-per ^ < gg > (mfe - < > 

2rr 4 



mj,mc — ^ — m^m^ 



1 , 9 1 n, fmbm^ — 1 (mc — + m? + 

' o"^o) ±r-^ + o^ 



1/ 



1 



+ 



.3 ^2 



— m^m^ 



_ ml (3mfe — mc) _ j 2{2mb — mc) mc 

(ml — Smlmc + Smtm'^ — mf — mbo"^ + mco'^) 1 ,^ 
+ ^ lt-2 \ ' (14) 

tyi ill ly I 

where r = — ml, r' = p'"^ — m^.. 

The QCD sum rules for the form factors f\{q'^) and f2{<f') are obtained by 
equating the phenomenological and QCD parts of the correlator and applying 
double Borel transformations with respect to the variables and p''^ {p^ — >■ 
Mf,p''^ M|) in order to suppress the contributions of higher states and 
continuum: 

/^(g2) _ (mb + mg) (m, + m^) ^^2^jMf^m%jM^^ 



ds r ds'p,{s,s',q^)e-'/''"e-''/^"^ 



(27r)2 7(m,+m,)2 

+5(Mi')5(M2)nr"~^"''|, (15) 

where i = 1,2 and B{Ml)B{M2) denotes the double Borel transformation 
operator. In Eq. (|T5l) . in order to subtract the contributions of the higher 
states and continuum, the quark-hadron duality assumption is used, i.e., it 
is assumed that 

pMgherstates^^^ s') = p^^^(s, s')e{s - So)e{s' - s'^) . (16) 
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In calculations, the following rule for double Borel transformations is also 
used: 

B — - > (-1)"^+^— ^p-rnl/Ml -m'i/Ml I q^n 

r-r'" ^ ' r(m)r(ri) (Mf)— i(M|)"-i ' ^ ^ 

3 Numerical analysis 

The sum rules expressions for the form factors fi{q^) and f2{(f) show that 
the condensates, leptonic decay constants of Bq and Dq mesons, contin- 
uum thresholds sq and s'q and Borel parameters and M| are the main 
input parameters. In further numerical analysis, we choose the value of 
the condensates at a fixed renormalization scale of about 1 GeV [23]: < 
dd >=< uu >= -(240 ± 10 MeVf, < ss >= (0.8 ± 0.2) < uu > and 
rriQ = 0.8 ± 0.2 GeV^. The experimental values for the mass of the mesons, 
TTiD, = 1968.2 ± 0.5 MeV, mD,{mr)±) = 1869.3 ± 0.4 MeV, mo^imDo) = 
1864.5 ± 0.4 MeV niB, = 5367.5 ± 1.8 MeV, mB^{mBo) = 5279.4 ± 0.5 MeV 
and mB^{mB±) = 5279.0 ± 0.5 MeV [23] are used. For the value of the lep- 
tonic decay constants and quark masses, we use the following values in two 
sets: In set 1, we use the results obtained from two-point QCD sum rules 
analysis: /s, = 209 ± 38 MeV m, fn. = 294 ± 27 MeV [3\, /b, = /b„ = 
140 ± 10 MeV [26] and fo^ = fo^ = 170 ± 20 MeV [26]. The quark masses 
are taken to be nicifi = m^) = 1.275 ± 0.015 GeV, m,(l GeV) ~ 142 MeV 
mail GeV) ~ 5 MeV, m„(l GeV) ~ 1.5 MeV and = (4.7±0.1) GeV 
\. In set 2, the recent experimental values fo^ = 274 ± 13 ± 7 MeV [27] , 
fD+ = 222.6 ± 16.7l|^ MeV [28], /b+ = 176tll?^ MeV p] and lattice 
prediction for /^^ = 206 ±10 MeV [2H] are used. For heavy quark masses 
TTic = 1.25 ± 0.09 GeV and mf, = 4.7 ± 0.07 GeV [23] and for light quark 
masses the values at the scale /i = IGeV (the same as set 1) are considered. 
The continuum threshold parameters Sq and Sq are also determined from 
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the two-point QCD sum rules: Sq = (35 ± 2) GeV^ ^ and s'^ = 6 GeV^ 
[3]. The Borel parameters Mf and M| are auxihary quantities and therefore 
the results of physical quantities should not depend on them. In QCD sum 
rules method, OPE is truncated at some finite order, leaving a residual de- 
pendence on the Borel parameters. For this reason, working regions for the 
Borel parameters should be chosen such that in these regions form factors are 
practically independent of them. The working regions for the Borel param- 
eters Mf and M| can be determined by requiring that, on the one side, the 
continuum contribution should be small, and on the other side, the contribu- 
tion of the operator with the highest dimension should be small. As a result 
of the above-mentioned requirements, the working regions are determined to 
be 10 GeV^ < Mf < 22 GeV^ and 4 GeV^ < M| < 10 GeV^. 

In order to estimate the decay width of Bg Dglu it is necessary to 
know the dependence of the form factors fi{q^) and f2{<f) in the whole 
physical region < < {msg — rriD^y. The dependencies of the form 
factors can be calculated from QCD sum rules (for details, see For 
extracting the dependencies of the form factors from QCD sum rules, we 
should consider a range of where the correlation function can reliably be 
calculated. For this purpose we have to stay approximately 1 GeV^ below 
the perturbative cut, i.e., up to = 10 GeV"^. In order to extend our results 
to the full physical region, we look for parametrization of the form factors 
in such a way that in the region < g^ < 10 GeV"^, this parametrization 
coincides with the sum rules prediction. The dependence of form factors 
/i(g^) and /2(g^) on g^ for set 1 are given in Figs. 3 and 4, respectively. 
Our numerical calculations show that the best parametrization of the form 
factors with respect to g^ are as follows: 
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l + aq + (3q^ + -ff + Xq'^ ' 
where q = q^ /m\^. The values of the parameters /i(0), a, /3, 7 and A for 
setl are given in the Table 1. 

Now, we are going to calculate the total decay width for these transitions. 
The differential decay width is as follows: 

,2\ 2 



dq^ 



1 



1927r3m^ 



q- - ml 



\{2q' + r4) [\h{q')\\2m\ + 2m\-q' 
+ 2{m\ - m\)Re[h{q')f;{q')] + \h{q'W 

^ ,,.2 , _2^ 



(r + mj) r 



\h{q')\\<-rn 



Do 



+ 2{mX - miyRe[Mq')f;{q')] + \f2{q')\V]} ■ (19) 



Next step is to calculate the value of the branching ratio for these decays. 
Taking into account the dependencies of the form factors and performing 
integration over in Eq. (|T9l) in the interval mf < q^ < {rriBq — rriDgY 





f(0) 


a 




7 


A 




0.24 


-1.57 


1.66 


-10.43 


19.06 




-0.13 


-1.69 


0.11 


1.50 


-4.65 




0.52 


-1.49 


0.02 


0.93 


-3.76 


h{Bu ^ D„(1864)£z/) 


-0.29 


-1.69 


0.21 


0.90 


-3.38 


f\{Bd ^ Dd{l8m)^u) 


0.52 


-1.49 


0.05 


0.77 


-3.47 


f2{Bd ^ Dd{l8m)^u) 


-0.29 


-1.69 


0.16 


1.13 


-3.70 



Table 1: Parameters appearing in the form factors of the Bq —>■ Dqiu 
decays in a four-parameter fit for Mf = 15 GeV"^, M| = 6 GeV^ and setl. 
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and using the total life-time tb, = 1.46 x 10 ^^s [33], r^^ = 1.64 x 10 ^^s, 
tb^ = 1.53 X I0~^^s [23J and | Kfe |= 0.0416±0.0006 [M], the following results 
of the branching ratios for set 1 are obtained. 

B{B, Dsiu) = (2.8 - 3.5) x 10"^ 
B{Bd Ddy) = (1.8 - 2.4) X 10~^ 
B(fi„ ^ DJz/) = (1.6 - 2.2) X 10"^ (20) 

The result for Bg DsCv shows that this transition can also be easily 
detected at LHC in the near future. The measurements of this channel and 
comparison of their results with that of the phenomenological methods like 
QCD sum rules could give useful information about the structure of the Ds 
meson. 

At the end of this section, we would like to compare the present work 
results of the form factors and their limits at heavy quark effective theory 
(HQET) (for details see |9J) for two sets with the predictions of the lattice 
QCD [T31 dn] at zero recoil limit for B Diu. For this aim, we introduce 
the notations used in [131 IIS] equivalent to Eq. (j3]) 



(21) 



< D I C7^,6 I B >= y/niBrnD 

where h+ and /i_ are the transition form factors and v and v' are the four 
velocities of the initial and final meson states. The relations between our 
form factors with the and h_ are given as: 

{niB + mD)h+ - {niB - mD)h_ 



fi 



J. _ {jub + mD)h_ - {rriB - mD)h+ 
^ 2^/m^m^ 

In order to perform the heavy quark mass limit, we define the multiplication 
of the V and v' as 

w = vv'= ^^^l-'i" . (23) 
2mBm£) 

12 



At zero recoil limit, w = 1 and from Eq. fl23l) it is correspond to ~ 
11 GeV^ which lies in the interval mf < < {uib ~ niDf'. Table 2 shows 
a comparison of the form factors and their HQET limits in the present work 
and the lattice QCD predictions at HQET and zero recoil {w = 1) limits 
in the present study notations. From this Table, it is clear that there is a 





/i 


/2 


Present study-set 1 


1.29 ±0.15 


-0.83 ±0.10 


Present study (HQET )-setl 


1.24 ±0.12 


-0.68 ±0.08 


Present study -set 2 


1.10 ±0.14 


-0.72 ±0.09 


Present study (HQET )-sct2 


1.06 ±0.10 


-0.58 ±0.06 


Lattice QCD (HQET )[13j 


1.19 ±0.01 


-0.68 ±0.05 


Lattice QCD (HQET )116J 


1.16 ±0.03 


-0.56 ±0.05 



Table 2: Comparison of the form factors in the present work, their HQET 
limits and lattice QCD predictions at HQET and zero recoil {w = 1) limits 
in the present study notations. 

good consistency among the models especially when we consider the errors. 
Moreover, a comparison of our results for the branching ratio of the Bg 
Dqiu with the predictions of the CQM model [22] and the experiment [23] 
are also given in Table 3. Considering the uncertainties and intervals, this 
Table also shows a good agreement among the phenomenological approaches 
and the experiment. Furthermore, this Table indicates that the value of 
the branching ratio increases both in the present work and the experiment 
by increasing the mass of the q quark. The intervals and uncertainties for 
values in the present study are related to the uncertainties in the values of 
the input parameters as well as different lepton types (e,yU, r). Our results 
for set 1 and set2 show that the value of the branching ratio is sensitive to the 
uncertainties in the value of the leptonic decay constants as well as the heavy 
quark masses. The existing uncertainties in light quark masses for q = u and 
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d cases don't change the results but for g = s case, we see a variation about 
3°/o in the value of the branching ratio. 

In conclusion, the semileptonic Bq Dqiu decays were investigated in 
QCD sum rules method. The g" dependencies of the transition form factors 
were evaluated. Using the expressions for the related form factors, the total 
decay width and the branching ratio for these decays have been estimated. 
The results enhance the possibility of observation of the Dgiu at 

LHC in the near future. Finally, the comparison of our results with that 
of the other phenomenological approach, lattice QCD and experiment was 
presented. 





(Bs ^ DJi^) 


(5, ^ 


(5„ DJu) 


Present study-set 1 


(2.8-3.5) X 10"" 


(1.8 - 2.4) X 10-" 


(1.6 - 2.2) X 10-" 


Present study-set 2 


(3.0-3.8) X 10"" 


(1.5 - 2.2) X 10-" 


(1.3 - 2.0) X 10-" 


CQM model 


(2.73 - 3.0) X 10-" 


(2.2 - 3.0) X 10-" 


(2.2 - 3.0) X 10-" 


Experiment 




(2.15 ±0.22) X 10-" 


(2.12 ±0.2) X 10-" 



Table 3: Comparison of the branching ratios for Bq — Dqiu 
decays in QCD sum rules approach, the CQM model [22] and the 

experiment j23j . 
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Figure 3: The dependence of /i on at Ml = 17 GeV^ , Ml = 6 GeV^ , 
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Figure 4: The dependence of /a on at Mf = 17 GeV^, M| = 6 GeV^, 
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